A refined invariant subspace method and 
applications to evolution equations 

g : Wen-Xiu Ma 

■ Department of Mathematics and Statistics, University of South Florida, 
Oh! 

<J ■ Tampa, FL 33620-5700, USA 

in 



> 
oo 

in 
in 

o 



Abstract: The invariant subspace method is refined to present more unity and more di- 
versity of exact solutions to evolution equations. The key idea is to take subspaces of solutions 
to linear ordinary differential equations as invariant subspaces that evolution equations admit. 
A two-component nonlinear system of dissipative equations was analyzed to shed light on the 
resulting theory, and two concrete examples are given to find invariant subspaces associated 
with 2nd-order and 3rd-order linear ordinary differential equations and their corresponding 
exact solutions with generalized separated variables. 
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1 Introduction 

Exact solutions to differential equations are significantly important in exploring the nature 
of motion. Based on the classification of elementary functions, there are only three kinds of 
explicit elementary exact solutions, classified as soliton-like, positon-like and complexiton-like 
solutions [1]- [I], to differential equations, both linear and nonlinear. Hirota direct method 
presents Wronskian and Pfaffian formulations of solutions to soliton equations, leading to soli- 
tons, positons and complexitons [5j El [7] . For general nonlinear partial differential equations, 
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symmetry related methods (see, e.g., [8]-[TT]) provide powerful approaches to their exact so- 
lutions. Group-invariant solutions stemming from symmetries play crucial roles in studying 
asymptotical behavior, blow up phenomena and fractal properties of motion, and they can be 
used to justify numerical schemes of solving partial differential equations [T2l [131 E]. 

The invariant subspace method, recently proposed in [151 IB] , is one of powerful approaches 
for constructing exact solutions to nonlinear evolution equations. Various invariant subspaces 
defined through linear ordinary differential equations have been presented for solving specific 
nonlinear evolution equations (see [T7]-[IH] and references therein). Indeed, the invariant sub- 
space method generates many interesting exact solutions to nonlinear evolution equations in 
mechanics and physics and a systematical solution procedure was given by Galaktionov and 
Svirshchevskii in their book [T7] . 

In particular, Galaktionov [16] utilized the invariant subspace method to generate exact 
solutions to nonlinear evolution equations with quadratic nonlinearities, and showed that exact 
positive solutions to the quasi-linear heat equations 

, _4 . _i r 

Ut = (u 3 u x ) x — au 3 4~ ou 3 + cu, 

where a,b,c G R are constants, can be constructed through invariant subspaces of functions 
of the polynomial or trigonometric form, admitted by the spatial differential operator. Actu- 
ally, evolution equations that admit invariant subspaces can be defined to be symmetries of 
given ordinary differential equations [T9J |20] . Interestingly, the iV-soliton solutions to soliton 
equations such as the KdV equation, the mKdV equation, the nonlinear Schrodinger equation 
and the sine-Gordon equation, derived by the Hirota bilinear method [3 [21], are all in a linear 
space of exponential functions under change of variables |17j . 

The invariant subspace method was also used to construct exact solutions to systems of 
nonlinear evolution equations. On the basis of the existence of invariant subspaces that sys- 
tems of linear ordinary differential equations define, Qu and Zhu [22] classified the systems of 
nonlinear parabolic equations of the form 

u t = [f(u, v)u x + p(u, v)v x ) x + r(u, v), 
v t = [g(u, v)u x + q(u, v)v x ] x + s(u, v). 

Zhu and Qu [23] presented an estimation of maximal dimensions of invariant subspaces for two- 
component systems of nonlinear evolution equations, and Shen, Qu, Jin and Ji [24J generalized 



2 



this estimation to multi-component systems of nonlinear evolution equations, together with 
some classifications of the considered systems of nonlinear parabolic equations and computation 
of the exact solutions derived from the corresponding invariant subspaces. 

In this paper, we would like to refine the invariant subspace method by taking invariant 
subspaces as subspaces of solution spaces to systems of linear ordinary differential equations. 
Note that a solution to an nth-order ordinary differential equation may not satisfy another or- 
dinary differential equation of order less than n. Our idea will generalize the invariant subspace 
method from the point of view of unity and diversity of invariant subspaces and exact solutions. 
A two-component nonlinear system of dissipative equations is analyzed carefully and a set of 
sufficient and necessary conditions is presented for the existence of invariant subspaces. Two 
concrete examples illustrate the effectiveness of the resulting refined theory in presenting exact 
solutions with generalized separated variables. 

2 Refining the invariant subspace method 

2.1 Scalar case 

Let us consider a scalar evolution equation 

u t = F[u], 

where u = u(x, t) is a function of x, t G R and F is a differential operator of order m: 

F[u] = F(x, t, «,«!,-•• , u m ), Ui = — { u, i > 0. 
Let n > 1 be a given natural number. Take n linearly independent functions 

f 1 (x),f 2 (x),--- ,f n (x), 
and form an n-dimensional linear space 

n 

W n = C{h{x), f 2 (x), ■■■ , f n (x)} = {J^Cifiix) | d = const., 1 < i < n}, (2.3) 

i=i 

i.e., the linear span of fi(x), / 2 (x), • • • , f n (x) over 1R or C. 

Definition 2.1. A finite- dimensional linear space W n is said to be invariant with respect to a 
differential operator F, if F[W n ] C W n , i.e., F[u] G W n , Wu G W n . 
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(2.1) 
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Suppose that W n is invariant with respect to a given differential operator F. Then there 
exist n functions F±, F 2 , ■ ■ • , F n such that 

n n 

F^QMx)} =Y,Fi(Ci,C 2 ,---,C n )f i (x) (2.4) 
i=i i=i 

for whatever constants C\, C 2 , ■ ■ ■ , C n . It follows that the evolution equation (12.11) possesses a 
solution of the form 

n 

u{x,t) = Y,<i> i (t)f i (x), (2.5) 
i=i 

if and only if 1; <p 2 , • • ■ ,<f> n satisfy a system of ordinary differential equations: 

-jr — Fii^U 02, • • ■ > 0n)> l<^<n. (2.6) 

We usually take an invariant subspace W 7 ™ as the space of solutions to a given nth-order 
linear ordinary differential equation: 

L[y] = + On-iCar)!/^- 1 ) + • ■ ■ + a (x)y = 0, y« = £fy, D = -^-, z>0, (2.7) 

where a ,ai,--- , a n _i are given continuous functions. The linearity of the above equation 
brings a good possibility to generate exact solutions to nonlinear evolution equations. 

The above approach for constructing exact solutions of the form (I2.5P is called the invariant 
subspace method [T5J [TB]. It is also called a generalized separation of variables [19J, whose 
resulting solutions of the form (12.51) we call solutions with generalized separated variables. 

To refine the invariant subspace method discussed above, let us consider a fc-dimensional 
subspace Wk of the n-dimensional linear space W n , and without loss of generality, we set 

k 

W k = £{A(x), f 2 (x), ■■■ , f k (x)} = QMx) | C t = const., 1 < i < k}, (2.8) 

i=i 

where k < n. The invariance condition Q Wk means that there exist k functions 

Fx, F 2 , ■ ■ ■ , F k such that 

k k 

F[J2CMx)] = ^F i (C 1 ,C 2 ,...,C h )f i (x) (2.9) 

1=1 i=l 

for whatever constants Ci,C 2 , - ■ ■ , C k . This way, beginning with a similar system of ordinary 
differential equations 

-|*=F i (V> 1 ,'0 2 ,..- ,Vfc), l<i< k, (2.10) 



we can engender a set of exact solutions to the evolution equation (12. ip : 

k 

u = Y^^ i {t)f i {x). (2.11) 

8=1 

We call this approach a refined invariant subspace method. One of its advantages is that when 
k < n, the invariance condition FfW^] C Wk requires much less conditions on the evolution 
equation (12. ip . 

The following results are helpful in searching for conditions to guarantee the existence of 
invariant subspaces that nonlinear evolution equations admit. 

Theorem 2.1. Let I = (a,b) be an open interval, xq £ I be an arbitrary point and m > be 
an integer. If a real function f defined on I satisfies 

/(z ) = 0, /'(z ) = 0, • ■ ■ , / (m - 1J (x ) = 0, f^(x ) = d, (2.12) 
where d is nonzero, then /,/',•"' > are linearly independent over I. 

Proof. If m = 0, the theorem is true, since based on (12.121) . / is not a zero function. In 
what follows, let m > 0. Suppose that the theorem is not true. That is to say, there are an 
integer < k < m and real constants q, < i < m — k, with c m -k 7^ such that 

c f(x) + Cl f'(x) + ■■■ + c m ^ k f m - k \x) = 0, x e I. (2.13) 

This is an ordinary differential equation with constant coefficients. Thus, the solution / is 
analytic in J, and so, using the conditions for derivatives in (I2.12p . we have 

f(x) = —Ax - x ) m + jr ^p^(x - x y. (2.14) 

t=m+l 

Now, the coefficient of (x — xo) k in the Taylor series expansion of the function on the left-hand 
side of (I2.13P is c m _kd/k\, which is not zero. This contradicts the linear dependence equation 
(I2.13p . Therefore, /,/',••• , f^ are linearly independent over /. □ 

Theorem 2.2. Let a,i(x), < % < n — 1, be real continuous functions on an open interval 
I = (a,b). Then there exits a solution y to the linear ordinary differential equation 

y^ + a n ^(x)y (n - 1) + --- + a (x)y = 

such that y,y', ■ ■ ■ ^y^ 1 " 1 ^ are linearly independent over I. 
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Proof. For a fixed xo G /, let us consider a Cauchy problem of the linear ordinary differential 
equation (I2.7P with the initial data: 

y(x ) = 0, y'(x ) = 0, • • • , y (n - 2) (xo) = 0, y^fa) = d, (2.15) 

where the real constant d is nonzero. The theory of ordinary differential equations tells that 
there exits a unique solution y defined on I to this Cauchy problem. It now follows from 
Theorem 12.11 with m = n — 1 that this solution y is the desired solution. □ 
Theorem 12.21 tells that for an nth-order linear ordinary differential equation, there always 
exists a solution y such that y,y', ■ ■ ■ , ?/ n_1 ) are linearly independent. Thus, in principle, for 
a given differential operator F, we can get sufficient and necessary conditions to guarantee 
the existence of an invariant subspace W kl by collecting all coefficients of linearly independent 
terms, generated from the linearly independent functions y,y',--- ,y ( - n ~ 1 \ in the invariance 
condition 

D n F[y] + an^D^Fiy] + ■■■ + a (x)F[y] = 0, y G W k , (2.16) 

and setting them to be zero. The process may be difficult and frustrating. However, from ana- 
lyzing different terms involving y,y', - ■ ■ , ?/( n_1 ) in f |2.16p . we can always get sufficient conditions 
for the existence of invariant subspaces that F admits. 

It should be interesting to note that W k may not possibly be generated by a fcth-order linear 
ordinary differential equation. An example is given as follows: 

W x = £{y}, y = e XlX + e X2X , A lj2 = consts., 

where y is a solution to the 2nd-order linear differential equation 

y" - (Ai + A 2 y + XxX 2 y = 0, 

but doesn't solve any lst-order linear differential equation when Ai ^ A 2 . Therefore, our 
refinement does make sense in generalizing the invariant subspace method. 

2.2 Multi-component case 

In what follows, we adopt the notations 

4 = u) = &U ^ X ; t \ l<i<q,j>l, (2.17) 
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such that the discussion can be easily extended to cases of multiple spatial variables. A system 
of evolution equations is assumed to take the form 

u t = F[u] = (*>], F 2 [u], F^[u}) T , u = {u\ u 2 ,---, ui) T , (2.18) 

where 

F*[u] = F\x, t, u\ • ■ ■ , u\ <), 1 < i < q, (2.19) 

are given sufficiently smooth functions in the indicated variables. Therefore, for each 1 < i < q, 
F % can be viewed differential operator of order rrij. 

Let Wk lt ...,k q denote a linear space W k \ x • • • x W^, with Wj:. being defined by 

Wl = £{fi(x), ■ ■ • , f ki (x)} = {J2 Cjfj(x) | C} = const., 1 < j < h}, 1 < * < q, (2.20) 

3=1 

where for each 1 < i < q, fl(x),--- , fl{x) are linearly independent. If the above vector 
differential operator F satisfies the invariance condition 

F[u] E W klr .., kq , Vw E W klr .. ikq , 

namely, 

F*[u] E Wl, Vu E W kw .., kq , l<t<q, (2.21) 

then the vector differential operator F (or the system of evolution equations (12.181) ) is said 
to admit an invariant subspace W klt ... tkq , or W klt ... tkq is said to be invariant under the given 
differential operator F. The above invariance condition (I2.2ip means that there exist functions 
Fj, 1 < j < h, 1 < i < q, such that 

fel kq kj 

^E^^)'--- 'E 6 '^)] =Y,*s(Clr~ ,Cl 1 ,...,Cl,...,Cl q )f j {x), (2.22) 

3=1 3=1 3=1 

where 1 < % < q. 

Now if a space W klr .. jkq is admitted by the vector differential operator F, then the system 
of evolution equations ( 12. 18ft possesses an exact solution of the form 

hi 

J2c;(t)f;(x), l<t<q, (2.23) 

3=1 



a 1 



if and only if the Cj(t)'s satisfy a system of ordinary differential equations: 



(2.24) 



The last step is that for each 1 < i < q, we take the space W].. = £{f{(x), ■ 
subspace of solutions to an njth-order linear ordinary differential equation: 



> fkS x )} as a 



^fo] = ir } + <-i(*)v< ( 



K-i) 



H h al(x)y- + ao(x)j/i = 0, 



(2.25) 



where rii > ki. The invariance conditions for the subspace Wk lt ...,k t — x • • • x with 
respect to F = (F 1 , • • • , F q ) T read 

£^[14] + ^(x)^- 1 ^] + ■ ■ ■ + a^F^u] =0, tiG W^,...,*,, 1 < 2 < g. (2.26) 

This set of equations is our starting point to construct exact solutions to systems of evolution 
equations by looking for their invariant subspaces. 

Note that the orders of linear ordinary differential equations defining invariant subspaces 
can not be arbitrary, and they are subject to the differential orders of the nonlinear operators 
F\ 1 < i < q. Once the maximal orders of the required linear ordinary differential equations 
are determined, we will be able to classify systems of evolution equations under consideration, 
and compute exact solutions from the associated invariant subspaces. 

The problem of maximal orders of linear ordinary differential equations defining invariant 
subspaces was firstly posed and solved for the scalar case in [IT] . For the scalar case, the 
maximal order of a linear ordinary differential equation defining an invariant subspace is not 
greater than 2m + 1, where m is the order of the differential operator F in (12. II) . For a q- 
component nonlinear really-coupled system defined by (12.181) and (12.191) with m\ > > • • • > 
m q > 0, the orders {ni, ■ • • ,n q } of linear ordinary differential equations defining invariant 
subspaces with rii > n 2 > • ■ • > n q > must satisfy [24J : 



That the system defined by (12 . 18[) and (12. 19[) is really-coupled means that for each pair 1 < i ^ 
j < q, there exists an integer < k < rrii such that is not a zero function. If F l , 1 < i < q, 
are all real, then the really-coupled condition can be concisely written as 




i=i 
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3 Invariant subspaces and exact solutions 

In this section, we analyze a (l+l)-dimensional nonlinear system of dissipative equations to 
illustrate how to generate invariant subspaces and the corresponding exact solutions. We con- 
sider the following nonlinear system of dissipative equations: 

u t = F = (u xx + aivv x ) x + a 2 v 2 , (3.1) 
v t = G = u xx + fan + f3 2 v, (3.2) 

where act, a 2 , Pi, (3 2 are constants, aci, a 2 are not simultaneously equal to zero, and we have used 
the traditional notation 

du dv d 2 u d 2 v 

V>x ~7{ ; V x — , U xx — -, V xx — -, • • • . 

OX ox OX A ox z 

Let us take an invariant subspace W 2j2 = W\ x W 2 defined by 

L x [y] = y" + a lV ' + a y = 0, L 2 [z] = z" + b x z' + b z = 0, (3.3) 

where a ,ai,b ,bi are constants to be determined. The corresponding invariance conditions 
read 

(D 2 F + aiDF + a oF)\ u&v ^ v( z W 2 = 0, (3.4) 
(D 2 G + b x DG + hG)\ ueW i v&W 2 = 0. (3.5) 

Substitute the expressions for F and G into the above equations, and replace u xx and v xx 
by — a\u x — ciqU and — b\V x — b^v a few times, respectively. Then we collect the coefficients of 
(v x ) 2 ,vv x and v 2 in the first simplified equation and the coefficients of u x and u in the second 
simplified equation, and set them to be zero, to obtain the sufficient conditions: 



v x ) 2 : 7aib 2 + a\a Q + 2a 2 — 4ai& — 3«iai6i = 0, 




(3.6) 


vv x : 12a 1 6 1 6 — aia bi — 4aiai6 + 2a 2 ai — ai&i 3 


— 2a 2 bi + aiai&i 2 = 0, 


(3.7) 


v 2 : 4ai6o 2 + ol 2 clq — a 1 6 1 2 6o + ai<iib bi — a\a bo ■ 


- 2a 2 b = 0, 


(3.8) 


u x : — ai 3 + 2 aofli — fii^i + a 2 b\ — ao&i + (i\b\ - 


- ai&o = 0, 


(3.9) 


u : — aoai 2 + ao 2 — /?iao + ao a i^i — a obo + /3i^o : 


= 0, 


(3.10) 



which guarantees the invariance conditions (13.41) and (13. 5p . We began with two second order 
differential equations, and so definitely there exist linearly dependent terms in (v x ) 2 ,vv x and 
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v 2 for whatever solution v, but u and u x could be linearly independent (see Theorem |2.2j) . 
Therefore, the conditions fl3.6p - fj3.8p are sufficient but not necessary to guarantee the first 
invariance condition (13.41) . but the conditions (13. 9p and (13.101) are both sufficient and necessary 
to guarantee the second invariance condition (13.51) . 

Under the first condition (I3.6p . the second and third conditions, (13.71) and (I3.8p . are equiv- 
alent to 

— aidibi 3 — 2 «i&o&i 2 + 3 «i6i 4 — a 2 a + a 2 a^b x = 0, (3-H) 
—3 «i6i 3 — a 2 ai + a\aibi 2 + 2 a^bo — 4 ai&o&i = 0. (3-12) 

This can be shown directly by using Maple and we will see later where they come from. 
Let us now assume A 2 = b 2 — 4&o > 0. Then 

W 2 = £{e x+x , e x ~ x }, (3.13) 

where 

A+ = - 



2 

Collecting the coefficients of three linearly independent terms e ( A ++ A -) x ; e ^+x anc [ e 2A_x j n ^ e 
first invariance condition (13. 4p and setting them to be zero gives rise to 

7i = 0, 72 ±^73 = 0, (3.14) 

respectively, where 

71 = aiao&i 2 — aidibi 3 — 2 a 2 aibi + «i6i 4 + 2 a 2 bi 2 + 2 a 2 a , (3.15) 

72 = —4 a 2 bo + a 2 a + 8 aibo 2 — 16 aibob\ 2 — 2a — la b + 4 a±bi 

+2 a 2 b 2 + 6 a\a\bf)b\ — 2 a\a\b\ — a 2 a\b\ + aiao&i 2 , (3.16) 

73 = —4 a\b\ — 2 a 2 bi + a 2 ai + 2 aia\bi 2 — 2 aiai&o — aiao^i + 8 ai&o&i- (3.17) 

The conditions in (13.141) equivalently lead to 

71 = 0,72 = 0,73 = 0. (3.18) 
These are sufficient and necessary conditions for guaranteeing the first invariance condition 

Under the first condition (13. 6p . 72 = and 73 = become the equations (13. lip and (13.121) . 
respectively. When the three equations, (I3.6p . 72 = and 73 = 0, hold, the condition 71 = 
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is automatically satisfied. Therefore, the conditions ( 13. 6p . (I3.7P and (13. 8p yield ( I3.18p . But 
conversely, it is not true. This is because the system (I3.18P has a solution a 2 = b = b x = 
with the other variables being arbitrary, but the left-hand side of (13. 6p under a 2 = bo = b x = 
is aido, not always zero. 

More specifically, under A 2 = b x 2 — 4&o > 0, let us take a smaller invariant subspace: 

W 2 ,i = W\ x W%, Wl = C{d lx }, (3.19) 

where fi = A + (or /i = A_). Then the invariance conditions 

(D 2 F + a\DF + a oF)\ u£W i )V( z W 2 = 0, (3.20) 

(D 2 G + b x DG + b G) \ ueW i jV&v a = 0, (3.21) 

only require 

72 + v^Ts = (or 7 2 - v^Ts = 0) (3.22) 

plus ( 13. 9 p and (I3.10p . Any of these two conditions in ( I3.22p is much weaker than the conditions 
in (I3.14p . i.e., ( I3.18p . Therefore, we can have a more general system of dissipative equations 
which still possesses exact solutions with generalized separated variables. 

In what follows, we give two concrete examples of getting exact solutions with generalized 
separated variables. 

Example 1: Let us consider a system 

9 1 

ut = (uxx + Oi X vv x ) x + (Sa^ibj - -&ibi 2 - -aia^v 2 , (3.23) 
v t = u xx + fcu + f3 2 v, (3.24) 

which admits an invariant subspace W 2 ,2 defined by 

L x [y] = y" + a iy ' + (a^ - b^y = 0, (3.25) 
L 2 [z] = z" + b x z' + { ai h - ~h 2 - ^ ai 2 )z = 0, (3.26) 

where a x can take any of the following three choices: 

, 4 6i 2 -A± y / 45 1 4 -20/3 1 6 1 2 +^7 

ai = b x or a x = — . (3.27) 

6oi 

We analyze the case of a x — b\, for which we have 

iH = (u xx + a x vv x ) x - 2a x a x 2 v 2 , (3.28) 
11 



v t = Uxx + 0iu + 2 v, (3.29) 

and 

L 1 [y]=y" + a 1 y' = 0, (3.30) 

L 2 [z] = z" + a lZ ' = 0. (3.31) 

From these two equations Li[y] = and l?^z\ = 0, we obtain an invariant subspace 

W 2 x Wf = £ {l, e" 01 *} x £ {1, e" aix } (3.32) 

that the system of (I3.28P and (I3.29P admits. It then follows that an exact solution takes the 
form 

u = d(t) + C 2 (t)e-^ x , v = D^t) + D 2 (t)e~ aix . (3.33) 

Substituting this solution into the system of (13.281) and (I3.29p . we get the following system of 
ordinary differential equations: 

C[ = -2a 1 a 1 2 D 1 2 , C 2 = -a^C 2 - 3a 1 a 1 2 D 1 D 2 , (3.34) 

D[ = 1 C 1 + 2 D U D' 2 = ai 2 C 2 + 0]C 2 + (3 2 D 2 . (3.35) 

Based on the refined theory, let us further focus on a smaller invariant subspace 

W 1A = Wl x W 2 = £{e- aix } x £{e~ aix }, (3.36) 

to present exact solutions. Solving the system of fj3.34j) and f!3.35|) with C\ — D\ — 0, we arrive 
at 

C 2 = ce~^\ D 2 = d<** - < a \ + M e -*\ (3.37) 

ai 6 + p 2 

and then according to (I3.33p . we obtain an exact solution to the system of (13.281) and fl 3 . 2 9 p : 

u = ce- aix - a ' h , v = de- aix+kt - c ( a \ + PJ e -a lX -a,H (3>38) 

ai 3 + P 2 

where c and d are arbitrary constants. 

Example 2: Let us finally consider a system of the following evolution equations: 

u t = (u x + aivv x ) x - a 2 a x v 2 - ?>a,\vv xx , (3.39) 

v t = u xx + a 2 u + 2 v, (3.40) 
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where a x ^ 0,ai,/?2 are arbitrary constants. This system admits an invariant subspace Wi t i 
defined by 

W X)X = Wl x Wl = £{cos(a lX )} x £{1 + sin(a 1 x)}. (3.41) 
These two basis solutions y = cos(aix) and z = 1 + sin(aix) satisfy 

L 1 [y]=y' / + a 1 2 y = 0, (3.42) 
L 2 [z] = z'" + a x 2 z' = 0, (3.43) 

but they cannot satisfy any lower-order linear ordinary differential equations with constant 
coefficients. 

Now, assuming a solution with the form 

u = C(t) cos(a x x), v = D(t)[l + sin(aix)], (3.44) 

and substituting back into the system of (13.391) and (13.401) . we find 

C\t) = - ai 2 C{t), D\t) = f3 2 D{t). (3.45) 

The general solution of this system yields the following exact solution to the system of (I3.39P 
and (ET20I1 : 

u = ce~ aiH cos(a 1 x), v = de p2t [l + &m{a x x)\, (3.46) 

where c and d are two arbitrary constants. 

This gives us a concrete example which uses the refined invariant subspace method to 
construct exact solutions to nonlinear systems of evolutions equations. 

4 Concluding remarks 

The invariant subspace method was refined by taking invariant subspaces as subspaces of so- 
lutions to linear ordinary differential equations. Our discussions were concentrated on how to 
identify sufficient and necessary conditions for the existence of invariant subspaces that non- 
linear evolution equations admit. Two concrete examples illustrated the effectiveness of the 
refined approach for exploring solution structures of systems of nonlinear differential equations, 
which are notoriously more difficult to solve than scalar ones. 

The invariant subspace method is also called a generalized separation of variables for nonlin- 
ear differential equations by Svirshchevskii [19] , presenting a kind of complexiton-like solutions 
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[2j [25] . It is interesting to see that the linear superposition principle takes on key role in con- 
structing exact solutions to either evolution equations [T7J or Hirota bilinear equations [26J. 
All related theories furnish linear combination solutions of functions with separated variables, 
which shows a sort of integrability of nonlinear differential equations [27]. 

Motivated by the multiple exp-function method [28], we can also characterize invariant 
subspaces through linear partial differential equations. This will lead to more diverse situations 
of solutions with generalized separated variables and open a much larger research area. 
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